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ABSTRACT 

We compute the acceleration on the Local Group using 11206 IRAS galaxies 
from the recently completed all-sky PSCz redshift survey. Measuring the ac- 
celeration vector in redshift space generates systematic uncertainties due to 
the redshift space distortions in the density field. We therefore assign galax- 
ies to their real space positions by adopting a non-parametric model for the 
velocity field that solely relies on the linear gravitational instability and lin- 
ear biasing hypotheses. Remaining systematic contributions to the measured 
acceleration vector are corrected for by using PSCz mock catalogues from 
N-body experiments. 

The resulting acceleration vector points ~ 15° away from the CMB dipolc 
apex, with a remarkable alignment between small and large scale contribu- 
tions. A considerable fraction (~ 65%) of the measured acceleration is gener- 
ated within 40 h~ x Mpc with a non-negligible contribution from scales between 
90 and 140 hr 1 Mpc after which the acceleration amplitude seems to have con- 
verged. The local group acceleration from PSCz appears to be consistent with 
the one determined from the IRAS 1.2 Jy galaxy catalogue once the different 
contributions from shot noise have been taken into account. The results are 
consistent with the gravitational instability hypothesis and do not indicate 
any strong deviations from the linear biasing relation on large scales. 

A maximum-likelihood analysis of the cumulative PSCz dipole is per- 
formed within a radius of 150 /i -1 Mpc in which we account for nonlinear 
effects, shot noise and finite sample size. The aim is to constrain the (3 = 
£l°- 6 /b parameter and the power spectrum of density fluctuations. We obtain 
(3 = OJOtoJ 5 at 1 a confidence level. 

The likelihood analysis is not very sensitive to the shape of the power spec- 
trum due to the rise in the amplitude of the dipole beyond 40 hT 1 Mpc and 
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due to the increase in shot noise on large scales. There is however a weak indi- 
cation that within the framework of CDM models the observed Local Group 
acceleration implies some excess power on large scales. 

Key words: Cosmology: theory - galaxies: clustering, - large-scale structure, 
large-scale dynamics. 
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1 INTRODUCTION 

If the Cosmic Microwave Background [CMB] defines a cosmological frame, then the dipole pattern observed in its 
temperature is a direct measure, via Doppler shift, of the Local Group [LG] velocity: v c = 627 ± 22 kms^ 1 towards 
(I, b) = (276° ±3°, 30° ±2°) as inferred from the 4-year-COBE data (Lineweaver et al. 1996). The best reason to prefer 
a velocity interpretation to more exotic scenarios (e.g. the Tolman-Bondi cosmological model of Paczynski & Piran 
1990, where the dipole moment of the CMB is of cosmological origin) is the remarkable alignment of the dipole with 
the gravitational acceleration vector measured at the LG location. This gravitational acceleration is inferred from the 
distribution of luminous objects in our cosmological neighbourhood and is expected to align with the dipole vector 
within the Gravitational Instability (Peebles 1980) and linear biasing framework in which the observed structures 
in the universe grew via gravitational instability from initial small fluctuations in the mass density field. In the 
linear regime, i.e. on scales where the mass density contrast, S p , is small, we also expect a direct proportionality 
between the amplitudes of the two vectors. The proportionality constant is a measure of the so-called /?(= 0.™ /b) 
parameter, where f2 m represents the density and b, the bias parameter, relates the mass density contrast, 5 P , to the 
fluctuations in the number density of the luminous objects, 5 g — b8 p . A comparison of the observed LG velocity with 
the gravitational acceleration is therefore a direct measure of /3. 

The many different estimates of the LG acceleration (see Strauss and Willick (1995) for a complete review) 
are partly contradictory. Most of the measurements agree in showing a small (< 30°) misalignment between the 
CMB and the LG dipoles. The amplitudes, however, are highly sensitive to the objects used as mass tracers. As 
reported by Strauss (1997), the LG acceleration as calculated using galaxies extracted from the IRAS 1.2 Jy redshift 
survey (Fisher et al. 1995) and from the Optical Redshift Survey [ORS] (Santiago et al. 1995, 1996) seems to receive 
little contribution from scales larger than 40 K~ x Mpc. Other analyses carried out using a deeper but sparser redshift 
survey of IRAS galaxies, known as the QDOT catalogue, have shown a non-negligible contribution from scales up to 
~ 100 h^ 1 Mpc (Rowan-Robinson et al. 1990). Clusters of galaxies, both optically and X-ray selected, can be used to 
probe very large scales because of their large luminosity. The LG gravitational acceleration generated by their spatial 
distribution also provides evidence of non-negligible contributions from depths up to ~ 150/i _1 Mpc (Scaramella, 
Vettolani & Zamorani 1991, Plionis & Valdarnini 1991, Plionis & Kolokotronis 1998). 

Such an apparent dichotomy may arise from the limited depth of the ORS and IRAS 1.2 Jy catalogues, from 
the sparse sampling in the QDOT and cluster catalogues, or from inhomogeneity and observational biases that may 
affect the different samples. 

In this paper, we use the recently completed PSCz survey of IRAS galaxies. This is a homogeneously selected all- 
sky catalogue which combines the depth of the QDOT with the dense sampling of the 1.2 Jy catalogue and therefore 
represents the ideal tool for investigating the large scale contribution to the LG gravitational acceleration. Our main 
concern is the combination of information on the growth of the dipole with other measurements of the local velocity 
fields (mainly from POTENT data) in a likelihood analysis and we therefore do not use the catalogue to its full 
extent. Rowan-Robinson et al. (1998) have done a careful analysis of the growth of the dipole on the largest scales 
using the same catalogue out to much larger distances. 

Even with an ideal flux-limited catalogue, however, other effects (shot noise, finite depth of the sample, use of 
the linear approximation, incomplete sky coverage, and redshift space distortions) cause the LG acceleration to differ 
from the observed LG velocity, making it impossible to determine f3 by direct comparison. In this work, we adopt an 
alternative two-step procedure in which: 

• the LG cumulative acceleration is computed from the data and its random and systematic errors are evaluated 
using PSC« mock catalogues extracted from N-body simulations. These catalogues are also used to statistically correct 
the measured PSCz dipole amplitude for systematic errors and to obtain a preliminary estimate of f3 from the data. 

• the resulting bias-free gravitational acceleration is then analysed using a likelihood approach similar to the one 
developed by Strauss et al. (1992) [S92 hereafter]. This will yield an estimate of /3 within the framework of a given 
cosmology and a way of evaluating errors. 
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The outline of this paper is as follows: In §2 we give a short theoretical introduction. §3 is devoted to computing 
the LG acceleration from the data and to estimating numerical biases and errors with the help of mock catalogues. 
§4 contains a likelihood analysis (theory and application) of the results, and conclusions follow in §5. 

2 THEORETICAL BACKGROUND 

In linear theory, the peculiar velocity, v, of a given galaxy is directly proportional to the gravitational acceleration, 
g, caused by the surrounding matter. The two vectors are therefore parallel and related by 



where / ~ fi^i . Hence, if the peculiar velocity is known, Q m can be determined by calculating the gravitational 
acceleration from the observed mass distribution, where Q m is the density of all matter in units of the critical density. 

During the past few years, particular attention has been paid to applying eqn. (|l]) to compute the acceleration 
acting on the Local Group of galaxies. The Local Group velocity in the CMB frame is accurately known from the 
CMB dipole anisotropy Q while the gravitational acceleration vector can be measured from redshift surveys which 
give an estimate of the three-dimensional galaxy density field. The vector g is related to the mass density field by 



where <5 p (r) is the density contrast at comoving coordinate r. Since we sample the density contrast of the galaxy 
rather than of the matter distribution, we need to postulate some relation between the two. In this work, we assume 
linear biasing, in which S g — bS p , such that we can shift from galaxy to matter density simply by dividing by the bias 
factor b. Eqn. (|l|) then becomes 



where f3 = f/b. 

In principle, /3 can easily be obtained by comparing |vy| and |v c |. In practice, however, the right hand side of 
eqn. (j^) is difficult to evaluate from redshift surveys because of several effects: 

• Finite Sample Size: The integral of equation |^ extends over all space. When using a redshift survey, however, this 
necessarily turns into a sum over a finite number of galaxies. The volume sampled extends to the maximum depth 
of the survey and excludes possible unsurveyed regions. Evaluating eqn. (^) in a finite survey volume results in a 
discrepancy between v r and v c . 

• Shot Noise: At large radii, the sampling of the galaxy distribution becomes more and more sparse, which leads 
to an increase in the so-called shot noise error. In this work, the shot noise contribution is evaluated following the 
S92 approach, which also accounts for the uncertainty in the galaxy mass function, i.e. the uncertainty associated 
with the assumption that all galaxies have the same mass. 

• Redshift Space Distortions: The surveys do not give full three-dimensional positions for the galaxies but rather 
angular positions and redshifts. Redshifts are related to distances via Hubble's law but only in the limit of very large 
distances. For nearby objects the peculiar velocities of the galaxies will add a significant contribution to the redshift 
that systematically distorts any density field calculated in redshift space (Kaiser 1987, Kaiser and Lahav 1988). We 
attempt to minimise redshift space distortions by reconstructing galaxy positions in real space from their redshifts. 

• Nonlinear Effects: Eqn. (|^) is not valid over the entire survey volume since linear theory breaks down in high 
density regions. Nonlinear contributions to peculiar velocities add incoherently to the vector, spoiling the linear 
relationship between velocity and acceleration vectors. 

* Note that both the CMB dipole and the gravitational acceleration caused by the surrounding matter are often referred to 
as 'dipoles'. In this paper we will obey this convention except when distinguishing between the 'CMB velocity' (v c ), which is 
inferred from the CMB dipole anisotropy, and the 'measured velocity' (v r ) inferred from the gravitational acceleration. 




(1) 




(2) 




(3) 
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All these effects have to be accounted for before a meaningful statement can be made about the value of 
calculated from a comparison of v,- and v c . In what follows, we will use a likelihood analysis that models the non- 
linear effects and the finite survey volume. The shot noise contribution will be calculated directly from the data 
while a self-consistent dynamical algorithm is adopted to model the errors introduced when minimising redshift space 
distortions. 

3 THE LOCAL GROUP ACCELERATION 

In this section, we present the dataset used, our methods for calculating the dipole and for evaluating the associated 
errors and finally the results of these procedures. By using simulated redshift catalogues, we hope to be able to correct 
the calculated dipole for all those errors that cannot be accounted for in the likelihood analysis. 

3.1 The Dataset 

Our dataset comes from the recently completed IRAS PSGz catalogue - a redshift survey of some 15,500 galaxies 
detected in the IRAS Point Source Catalogue. The PSOz catalogue contains almost every galaxy in the IRAS PSC 
with 60 (j,m flux (/6o) larger than 0.6 Jy for which the redshift has been measured. Only sources from PSC with 
/eo > O.5/25 (to exclude stars) and with /100 < 4/eo (to exclude most of the galactic cirrus) have been included in 
the catalogue. Our subsample contains 11206 PSCz objects within 200 h~ x Mpc , with PSC fluxes at 60 /im > 0.6 
Jansky, and with positive galaxy identifications. Regions not surveyed by IRAS (two thin strips in ecliptic longitude 
and the area near the galactic plane defined by a V-band extinction of > 1.5 mag.) are excluded from the catalogue, 
which therefore covers ~ 84% of the sky. 

We also consider a catalogue of galaxies extracted from the 1.2 Jy redshift survey (Fisher et al. 1995) which 
contains all the objects with /60 > 1.2 Jy, also within 20,000 kms -1 . The resulting 1.2 Jy subsample has ~ 88% sky 
coverage and contains 4626 galaxies. 

3.2 Minimising Redshift Space Distortions 

As outlined by Kaiser (1987) and Kaiser and Lahav (1988), redshift space distortions can modify the amplitude of 
the galaxy dipole, especially when it is measured from a flux-limited sample. To minimise this effect, we compute the 
galaxies' peculiar velocities and real space positions iteratively starting from the galaxy distribution as seen in the 
redshift space. A detailed methodological description can be found in Branchini et al. (1998), here we just summarise 
the main concepts. The method described is very similar to that developed by Yahil et al. (1991). 
We assume linear GI and LB and look for an iterative solution to the system of equations 



where V; is the peculiar velocity of the generic object i, v c is the LG velocity (both velocities are measured in the 
CMB frame), f is the unit vector along the radial direction. Distances and velocities are expressed in kms -1 . The 
redshift z; is measured in the LG frame. First, we make a guess for r (which in the 0-th iteration coincides with the 
galaxy redshift), then compute v(r) from eqn. (^), correct to a new r by eqn. (|E|) and continue the process until the 
calculated velocities Vi converge. This is typically achieved within 8 iterations. 

There are several additions to this, designed to minimise the uncertainties inherent in the procedure. When 
computing the gravitational acceleration from eqn. (^) we smooth the density field using a Top Hat window and a 
smoothing length equal to the average separation of objects at that distance. The minimum smoothing length is 5 
h~ l Mpc to ensure sufficient smoothing even at distances where the interparticle separation is small. This procedure 




(4) 



and 



CZi - f • (v, - v c ), 



(5) 
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mitigates nonlinear effects especially in high density environments but also impairs our ability to predict the small- 
scale features in velocity and density fields. 

The free parameter (3 has to be determined a posteriori by comparing observations with model predictions. We 
adiabatically increase f3 from 0.1 to 1.0 in 10 iterations, thereby slowly increasing the effect of gravity to improve 
convergence of the velocity field (S92). The resulting peculiar velocities scale to first order in (3 which, as we will 
discuss in §3.4, causes the galaxy dipole to be almost independent of the /3 value used in the reconstruction. 

The iterative procedure uses all the galaxies within 200 Mpc while the mass distribution beyond is considered 
to be homogeneous, i.e. we neglect its gravitational effect. This assumption is justified as long as the reconstruction 
is performed in the LG frame since relative peculiar velocities are not affected by external dipoles. Higher moments 
of external force fields only affect the reconstruction reliability in the outer regions. To avoid these edge effects we 
calculate peculiar velocities and real space positions only for galaxies out to 170 ft -1 Mpc even though material out to 
200 h^ 1 Mpc is considered to be gravitating. The error analysis in §3.4 indicates that the combined effect of shot noise 
and reconstruction uncertainties at that distance have already increased to over 30 % of the measured acceleration, so 
to maximise the reliability of the data, we perform the likelihood analysis only on data measured out to 150 h^ 1 Mpc . 

Spurious gravitational acceleration may also arise if the galaxy sample does not have full sky coverage. To alleviate 
this effect, we fill the unobserved regions with synthetic objects using two distinct procedures: At high (\b\ > 8°) 
galactic latitude we fill the masked areas with a uniform distribution of simulated objects having the same average 
number density and selection function as the PSCz galaxies. At lower galactic latitudes (\b\ < 8°, the so called 'zone 
of avoidance') the masked area is filled with galaxies which are cloned from the PSCz galaxy distribution in two 
strips above and below the zone of avoidance (Yahil et al. 1991), also retaining PSCz galaxies with |&| < 8° at their 
observed position. 

Since we are using a flux-limited catalogue, the number density of galaxies within the survey volume decreases 
with the distance. In the reconstruction we account for this effect by assigning a weight to each galaxy that is 
proportional to the inverse of the selection function at the position of the galaxy. Since the positions of the galaxies 
vary during the iterations, we recompute the selection function at each step from the observed galaxy fluxes using a 
parametric maximum likelihood technique, and the galaxy weights are updated accordingly. For the selection function 
we assume the analytic form proposed Yahil et al. (1991). 

Most of the results presented here refer to this iterative reconstruction procedure. However, two other methods, 
fully described and tested in Branchini et al. (1998), have been implemented to check the robustness of the results. 
The first is a grid-based version of the above iterative method while the second is based on the Nusser and Davis 
approach (1994) that uses the first order Zel'dovich approximation to correct for redshift space distortions without 
any need for iterations. In what follows we will refer to these three methods as Ml, M2 and M3. As we will show in 
§3.5, the PSCz dipoles reconstructed with the three methods are fully consistent with each other. 

3.3 PSCz Mock Catalogues 

Systematic uncertainties affecting the measured PSCz dipole have been evaluated using a series of PSCz mock 
catalogues. These were extracted from N-body simulations of CDM universes performed by Cole et al. (1998) using 
a 345.6 Mpc computational box. For the present analysis we have considered two of their cosmological models. 
They are: 

• an Q. m = 1 CDM model with T = 0.25 

• an Q A = 0.7, fi m = 0.3 CDM model 

These two are flat CDM cosmologies normalised to reproduce the local abundance of clusters as explained in Eke, 
Cole, and Frenk (1996): 

a 8 = 0.52n-°' 46+0 ' lonm if fl A = 0, and 

a s = 0.52n^°' 52+013nm if Q A = 1 - H m . 



© 0000 RAS, MNRAS 000, 000-000 



Local Group Acceleration 7 



The velocity fields of the N-body models have been smoothed with a top hat function of 1.5 h^ 1 Mpcto obtain 
pairwise peculiar velocities of ~ 250 kms _1 at 1 h^ 1 Mpc , a value close to that observed in the galaxy distribution 
(Guzzo et al. 1997, Strauss, Ostriker and Cen 1998). For each of the two models we produce 10 mock catalogues: We 
impose b = 1, a value compatible with the one commonly accepted for IRAS galaxies, in order to identify galaxies 
with particles in the simulations. The velocity field in each catalogue mirrors the conditions in the local universe: 
the central LG-like observer has a peculiar velocity of v c = 625 ± 25 kms -1 , the LG velocity residuals ((vi — w c ) 2 ) ' 5 
within 5/i _1 Mpc (sometimes called "shear") are smaller than 200 kms~ 1 and the fractional overdensity within the 
same region is in the range —0.2 < 5 < 1.0 (Brown and Peebles 1987). Each catalogue has a depth of 170/i _1 Mpc 
(since our velocity analysis of the PSCz will only extend out to that distance) and the coordinate system is rotated 
such that the observer's velocity vector points towards the observed CMB apex. Since the number density in the 
simulations (~ 0.039 h 3 Mpc -3 ) is smaller than the one given by the PSCz selection function within 10.9 h^ 1 Mpc, 
we are forced to semi-volume limit our catalogues at that radius. At greater distances we use a Monte Carlo rejection 
technique to force the simulated galaxy population to obey the N(z) distribution of PSCz galaxies. A random flux 
consistent with the selection function is then attributed to each galaxy. Finally, we mask the same areas as in the 
PSCz catalogue to mimic the incomplete sky coverage. A similar set of catalogues has been produced to mimic the 
galaxy distribution in the IRAS 1.2 Jy subcatalogue (§3.1). Given the smaller number density of 1.2 Jy galaxies 
(Fisher et al. 1995), the 1.2 Jy mock catalogues are semi-volume limited to a smaller radius of 7.8/i _1 Mpc. 

3.4 Error Estimates Using Mock Catalogues 

Several different effects prevent the reconstruction procedure from recovering the exact real space galaxy distribution 
and a less than perfect reconstruction generates both systematic and random errors in the dipole. S92 showed that 
the most serious systematic bias originates from the so called 'rocket effect' (Kaiser 1987, Kaiser and Lahav 1988) 
which represents the spurious acceleration measured from a magnitude-limited sample of galaxies by an observer 
with some peculiar velocity unconnected to the gravitational acceleration (e.g. an observer in a rocket). In our case 
the errors in the computed LG velocity mimic the non-gravitational acceleration. The characteristic signature of this 
effect is a monotonically increasing spurious contribution to the cumulative LG acceleration. S92 correct for this 
effect by modifying their reconstruction procedure. Here we adopt a different strategy which uses mock catalogues 
to correct not only for the rocket effect but also for further possible systematic effects that arise from filling the 
unsurveyed areas with synthetic objects. We consider a generic PSCz mock catalogue and measure the gravitational 
acceleration at the observer's central position from the true, all-sky galaxy distribution in real space. We then measure 
the gravity vector from the same catalogue, after applying the filling procedure and after iteratively moving galaxies 
to their reconstructed real space positions. The two acceleration vectors are similarly affected by shot noise, finite 
sample size and nonlinear effects and thus any discrepancy can only be ascribed to the intrinsic uncertainties in 
the reconstruction, the filling techniques and to the rocket effect. For each of the mock catalogues we compute the 
vector difference between the two cumulative accelerations at different radii, thereby obtaining an estimate for the 
systematic dipole error. We then average the difference over the various catalogues of both cosmologies to account for 
our ignorance of the underlying, true cosmological model. The mean discrepancy represents the cumulative effect of 
the various systematic errors while the dispersion around the mean quantifies the random errors. A similar exercise 
has been performed using mock IRAS 1.2 Jy catalogues. Results of this procedure are displayed in figure 1. The points 
represent the systematic error on the cumulative dipole measured at a given radius whereas the errorbars display the 
dispersion in that error. Since we are showing a cumulative quantity, the errors on each shell are correlated, both 
in amplitude and direction. Note that Vtrue is the dipole velocity calculated using real space positions of galaxies, 
whereas v rec is computed using the redshift space positions as described in §3.2. Both quantities are normalised to 
P — 1. Filled dots refer to PSCz mocks while small triangles refer an analogous analysis performed using the simulated 
IRAS 1.2 Jy catalogues. 

It is encouraging to note that there is virtually no systematic error up to a distance of 150 h^ 1 Mpc . The marginal 
increase beyond that is probably the signature of the Kaiser effect. In what follows, the average systematic offset will 
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Figure 1. Error estimates from the mock catalogues. Black dots show the discrepancy between the true and the "observed" 
PSCz dipole, averaged over different mock catalogues at different radii. Error bars represent the scatter around the mean and 
quantify the amplitude of the random errors. Small triangles show similar results for the 1.2 Jy mock dipoles. 



be subtracted out of the PSCz and IRAS 1.2 Jy dipoles to statistically correct for rocket-like effects. The importance 
of the different error sources seems to vary with the distance. In the inner regions nonlinear effects and intrinsic error 
in the reconstruction seem to dominate. Uncertainties in the filling procedures and Kaiser effect becomes dominant 
on larger scales. The former accounting for ~ 40 % error budget at a distance of 100 h^ 1 Mpc , both in amplitude and 
direction. The typical random error on the differential dipole, a, amounts to ~ 15 kms _1 per shell of 10 h^ 1 Mpc. 
This error will be accounted for in the likelihood analysis of §4 in the same way as the shot noise. 

The errors showed in figure 1 turned out to be Gaussian distributed around their mean (i.e. the black dots) which 
indicates that our error estimate, albeit based upon only 20 PSCz mock catalogs, are realistic. As outlined in §3.2, 
the reconstructed dipole scales with {3 linearly only to a first approximation. The actual result marginally depends 
on the value of j3 assumed in the reconstruction itself. To evaluate the amplitude of this effect we have performed 9 
PSCz reconstructions varying f3 in the range [0.4,1.2], to include the range allowed by observations (see Dekel 1997). 
Figure 2 shows the average cumulative dipole from the 9 iterations. The upper and lower dashed lines display the 
dipoles for the f3 = 1.2 and /3 = 0.4, respectively. At the radius of 150 /i -1 Mpc, the percentage difference is 11 %, 
notably smaller than the 32 % caused by the cumulative effect of random errors and shot noise. Also, note that the 
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Figure 2. Real Space PSCz cumulative dipole reconstructed with the Ml method using different a priori values of /3 between 
0.4 and 1.2. Solid line is the average of all results, dashed lines represent the dipoles calculated with the highest and lowest 
values for (3. 

shape of the dipole remains constant at large radii which. This corroborates our previous conclusion that systematic 
errors induced by the Kaiser rocket effect are far less prominent than in the S92 analysis. 

Uncertainties deriving from the nonlinear f3 scaling will therefore be neglected in the rest of the work. 

3.5 The PSCz Dipole 

We are now ready to compute the PSCz dipole corrected for systematic biases. Figure 3 shows the magnitude of the 
cumulative PSCz acceleration up to the limiting depth of 170 /i _1 Mpc, computed after the reconstruction procedure 
but with no correction for systematic effects. The amplitude is normalised to (3 = 1.0. The real space dipole computed 
using Ml is shown by filled circles while the dipoles from M2 and M3 reconstructions are displayed with squares and 
open circles, respectively. Error bars representing the same random errors as in figure 1 are displayed for Ml alone to 
avoid overcrowding. The cumulative dipoles of the three different reconstructions appear to be remarkably similar, 
despite the very different nature of the techniques used. The same considerations also apply for the dipoles' directions 
and speak for the reliability of our dipole estimate. 

Since the three methods produce equivalent results let us, from now on, concentrate on the Ml dipole only. 

Table 1 shows the amplitude of the uncorrected cumulative dipole and the shot noise contribution, both normal- 
ized to (3 = 1, along with the dipole direction along I and b and the cumulative misalignment angle of the reconstructed 
PSCz dipole from the CMB dipole apex The latter quantity is shown in the left panel of figure 4. The misalignment 
is remarkably small already beyond 40ft~ 1 Mpc, i.e. when large overdensities (Great Attractor and Perseus Pisces 
supercluster) enter the sampled volume. This indicates that contributions to the LG accelerations are coherent over a 
very large range of scales (Basilakos and Plionis 1998). The right hand panel displays the direction of the cumulative 
gravity vector measured within 150/i _1 Mpc while the error bars quantify the random uncertainties derived from the 
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Distance (h _1 Mpc) 



Figure 3. PSCz cumulative acceleration: the amplitude. Black circles represent the amplitude, normalised to (i = 1, of the 
cumulative PSCz dipolc obtained when correcting for redshift space distortions using method 1 (squares = method 2, open 
circles = method 3). Error bars are the the same as in figure 1. The dashed line on the bottom represents the cumulative shot 
noise error. No correction for systematic uncertainties has been applied. 

mocks. The continuous curves are loci of constant misalignment in the (I, b) plane. The asymptotic misalignment 
angle between LG and CMB dipole 56 — 15° is remarkably small. 

Both PSCz and 1.2 Jy samples have been extracted from the same parent PSC IRAS catalogue and have a 
similar sky coverage. The main difference between the two is the limiting flux which causes the radial selection to be 
more severe for the IRAS 1.2 Jy galaxies. The dipoles from the two samples should therefore exhibit similar features 
apart from different shot noise errors. As shown in figure 5, we indeed find a remarkable agreement. The left panel 
shows the cumulative amplitudes of the two dipoles, both normalised to f3 = 1. Filled dots and small triangles refer 
to the PSCz and 1.2 Jy samples, respectively. Error bars represent random errors for the PSCz dipole alone. Unlike 
in figure 4, dipoles in figure 5 have been statistically corrected for the systematic errors (as calculated in §3.4), which 
mainly arise from the rocket effect and the uncertainties in the filling procedure. The two curves on the bottom 
display the shot noise errors for the PSCz (continuous line) and 1.2 Jy (dashed line). Both corrected dipoles show 
the same features: a sharp rise followed by a small decrease around 50 h^ 1 Mpc as a result of the competing pull of 
the Great Attractor and Perseus Pisces regions. There is a further slow rise in the PSCz dipole of ~ 35% in the range 
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Figure 4. PSCz cumulative acceleration: the direction. The plot on the left shows the cumulative misalignment between PSCz 
and CMB dipole vectors. On the right the black dot represents the direction on the sky of the cumulative PSC dipole as 
measured at 150 ft -1 Mpc . Errorbars quantify 1— a random errors from the mocks. The starred symbol at the centre shows the 
CMB dipole directions. Contours are drawn at constant misalignment angles. 

between [60-140] h~ l Mpc. Both the PSCz and the 1.2 Jy show significant evidence for gravitational contributions 
beyond 40 /t -1 Mpc. However, unlike the amplitude of the 1.2 Jy dipole, that of the PSCz seems to be constant 
beyond ~ 140 h~ l Mpc. This discrepancy can be accounted for once the different shot noise contributions at large 
radii are taken into account. The right hand side of figure 5 shows how large the dipole's misalignment becomes when 
using the 1.2 Jy catalogue. The 1.2 Jy gravitational acceleration within 150 ft -1 Mpc points ~ 25° away from the 
CMB dipole, which reduces to a mere ~ 15° when the PSCz sample is considered (see also Kolokotronis el al. (1996) 
who predicted the misalignment of the PSCz dipole to lie between 10 and 17 degrees). Error bars represent random 
errors - no correction for systematic effects is needed as they mainly affect the dipole's amplitude. 

A further proof of the robustness of our result comes from the comparison with preliminary results from the 
analysis of the PSCz dipole by Rowan- Robinson el al. (1998), which used a different methodology and a much deeper 
PSCz subsample to address the issue of its convergence. The two dipole estimates agree well within 1—cr errors. Also 
note that the above results are fully consistent the recent re-analysis of the 1.2 Jy and QDOT catalogues by Basilakos 
and Plionis (1998). 
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Table 1. Cumulative PSCz Dipolc. Column 1: Distance from the LG in ; Column 2: Cumulative amplitude in (/3 = 1); Column 
3: Random errors in (/3 = 1); Column 4: Cumulative direction along I ; Column 5: Cumulative direction along b ; Column 6: 
Cumulative misalignment w.r.t. the direction of the CMB dipole. 



3.6 A Numerical Estimate of (3 

In this section, we will evaluate /3 in a purely numerical way to compare to the likelihood analysis. The idea is similar 
to that of §3.4, i.e. we use the mock catalogues to evaluate the errors in the reconstructed dipole. Here mock PSCz 
catalogues are used to statistically correct for the discrepancy between the measured acceleration vector and the 
true peculiar velocity of an observer. In this exercise we are assuming that the two cosmological models assumed 
to generate the mock catalogues are realistic. Indeed, the two cosmologies explored produce large scale motions 
compatible with the observational constraints (Jenkins et al. 1998). Moreover, for the cosmological models we have 
considered, the contribution to the dipole from the wavelength longer than the size of the N-body computational 
box is small compared to the others errors. Even more so in the inner regions, where a very significant contribution 
to the LG acceleration comes from. Therefore we do not account for this error source in the present analysis. 

For each mock catalogue we compute, at different radii, the ratio between the reconstructed dipole at the 
observer's position and its true N-body velocity. The average of the ratios calculated from the various mock catalogues 
provides us with a multiplicative factor which relates the dipole reconstructed by our method to the real LG velocity. 
We calculate this for all three Cartesian components and find that the multiplicative factors are sufficiently similar 
to warrant averaging. 

Note that, unlike in §3.4, the discrepancy between reconstructed dipole and real velocity is caused by the effects 
of non-linearity, shot noise, finite volume and residual redshift space distortions i.e. this now represents the cumulative 
effect of all sources of error. Note also that the multiplicative factor and its dispersion do not depend on /3 in the first 
approximation when it has been averaged over all catalogues. The reconstructed dipoles have been scaled to what is 
in each case known to be the correct /3 of the simulation so that the ratio between the reconstructed and the real 
dipole velocity only depends on the various errors and not on cosmology. 

We then use this multiplicative factor to recover the LG peculiar velocity from the reconstructed cumulative 
PSCz dipole to be compared with CMB dipole, v c , to evaluate /3. 

Figure 6 shows the result of this procedure. Surprisingly, /3 is remarkably stable around 0.7 even at low radii, 
where the volume includes little of the mass distribution that really causes the total LG acceleration. However, the 
errorbars on those scales are correspondingly larger. The errors, which represent the dispersion in the distribution 
of cumulative factors obtained from the mocks, are asymmetric, which indicates that it will be easier to set a lower 
rather than an upper limit on /3. 

The dashed and two dotted lines show the result of the likelihood analysis performed in §4.3 (figure 7a and 7b). 
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Figure 5. PSCz vs. 1.2 Jy dipoles. The amplitudes of the PSCz and the 1.2 Jy dipoles, both corrected for systematic errors, 
are shown in the left plot using filled dots and triangles, respectively. Error bars represent 1— a uncertainties. The lines on the 
bottom indicate the cumulative shot noise error for the PSCz dipole (continuous line) and for the 1.2 Jy one (dashed line). On 
the right plot the same symbols indicate the directions of the two dipoles at 150 h~ l Mpc . 

The dashed line represents the most likely /3-value, the dotted lines the 1 — a range. It is encouraging that the 
numerical analysis is entirely consistent with the likelihood analysis. 



4 LIKELIHOOD ANALYSIS 
4.1 Velocity Statistics 

Our aim is to assess the likelihood of a given cosmological model using all the information we have on velocities in 
the local universe, i.e. the measured LG velocity, the CMB dipole, the bulk velocity, and the LG velocity residual. To 
this end, we need to find a mathematical expression for the probability distribution functions (joint or constrained) 
of the velocities available and the first step is to find expressions for these velocities. In the present analysis, which is 
based on the work of Juszkievicz, Vittorio and Wise (1990) and S92, all velocities are regarded as Gaussian random 
vectors. 

In the limit of a finite survey volume, the peculiar velocity of the Local Group becomes 
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Figure 6. Robustness Test. Filled dots represent the /3 value obtained when comparing the LG velocity deduced from the CMB 
dipole to the one numerically deduced from the PSCz acceleration. Errorbars represent 1— a uncertainties estimated from the 
mock PSCz catalogs. The dotted and the dashed lines represent the 1— a level range for the value of beta from the likelihood 
analysis (J3 = 0.7jl o ;|). 

Vi = ^ / dWiM^p^r/) (6) 

where the window function Wi(r) specifies the finite survey volume. It is often more convenient to work in Fourier 
space in which case eqn. (^) becomes 

v, = H o p J d 3 k^(fc)^i^«5 s (k), (7) 

where k is the wavenumber vector, Wi(k) — k J drWi(r)ji(kr), is the window function in Fourier space and ji(kr) is 

the spherical Bessel function of order 1. This expression represents the linear velocity of the LG that is expected to 

be measured through a given window function in a universe characterised by the density field <5 9 . To make use of the 

information stored in the growth of the cumulative LG acceleration vector with radius of the sampled volume, we will 

find it more convenient to measure the differential LG velocity, v r ,j generated within a series of N non-overlapping 

top-hat window functions between radii fj and r^+i, 

w = Mkn) _ sin(fcr i+1 ) 
kri kr i+ i 
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(a) a B = 0.52, shot noise only 



(b) as (a) with a added 
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(f) as (c) with <r a = 0.87 




Figure 7. Probability of Amplitude and Misalignment. Likelihood contours for amplitude and misalignment of the PSCz dipolc 
in a specified world model. Fig 7a refers to the case of (F = 0.5,/3 = 1,(78 = 0.52) universe while contours refer to C(v r | ~vcmb)- 
Fig 7b show the effect of introducing the reconstruction errors, a, and Fig 7c also includes the decoherence effects. Fig 7d and 
c show the effect of introducing either the shear or the bulk flow constraints, respectively. Figure f show the same case as figure 
7d for a (r = 0.195, (3 = 0.55 and a 8 = 0.87) model. 



with i = 2, ...N. In the innermost region (i — 1) we smooth with a top hat filter on scales smaller than r s = 5 h x Mpc 
to eliminate strong nonlinear effects. Therefore, the window function for the innermost region is of the form 



kr s 



(9) 



The total measured LG velocity v r can be obtained when N = 1 and n = R max is the sample's depth. 

The other velocities modelled in this analysis are equally characterised by different window functions. The iV+lth 
velocity vector considered is the true LG velocity measured in the CMB frame, v c . Its window function extends over 
all space except for a small-scale cut-off to account for the finite size of the Local Group. In fc-space its form is 
sin(fcrLG) 



WcMB = 



(10) 



(kr L c) 
where tlg — 1 h~ 1 Mpc. 

The velocity field within the neighbourhood of the Local Group is remarkably quiet (c/.Peebles 1998). The 
observed LG velocity residuals v s , i.e. the difference between the Local Group velocity and the average velocity 
within a sphere of 5ft _1 Mpc, is less than 200 kms" . This can be used as an effective constraint in a likelihood 
analysis of the LG dipole field and we take the amplitude of the velocity residuals as the N + 2th random variable. 
The window function for the LG velocity residuals within r a — 5 /i _1 Mpc is 
sin(fcr a ) — kr a cos(fcr a ) 



W s {k) = 1-3- 



{kr a 



(11) 



The final N + 3th random vector we consider is the bulk velocity v;, measured within = 30/i _1 Mpc from the 
POTENT analysis of the Mark III catalogue of galaxy peculiar velocities (Dekel 1997). Since we consider POTENT- 
processed data the measured bulk velocity is just the average peculiar velocity of a spherical region of 30/i _1 Mpc. 
This allows us to express the relative window function for the bulk flow as 
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~ , Q sin(fcr h ) - kr b cos(kr b ) 

Wb ^ = 3 (kny* ■ (12) 

Since we model all the N + 3 velocities as Gaussian random vectors, any probability distribution function will 
only depend on the covariance matrix: 

M im = i < v; ■ v m > 



H 2 o p 2 
6tt 2 



dkP{k)W l {k)W m {k), (13) 



where P(k) = (8(k)8(k)) is the power spectrum of density fluctuations and Wi(k), W m (k) are any of the window 
functions that we have previously introduced. In this paper, we restrict our analysis to the family of CDM power 
spectra for which we use the expression of Davis et al. (1985), characterised by the shape factor T — Q, m h: 

P(k) = ^exp[-l(% 3 ) 2 ] 

v 1 [i + i.7fc/r + 9(fc/r)-V2 + (fc/r)2]2 ^ > 

where g s = 100 km s -1 and the only function of the exponential term is to improve numerical stability, h is the Hubble 
constant in units of 100 kms -1 Mpc -1 and B is the normalisation constant. As anticipated in §3.3, we normalise the 
power spectrum to as according to the Eke, Cole and Frenk (1996) prescription. Thus, we can completely characterise 
the CDM background cosmology by specifying /?, F and Qa- 



4.2 Probability Distributions 

The Gaussianity of the velocity field guarantees that the joint distribution function of all our N+3 model velocities 
is also a multivariate Gaussian: 

f(v r (i),i= l..JV,v c ,v s ,v 6 ) = (27r)- 3(JV+3)/2 (detM)- 3/2 exp(-ivi • v m (M- 1 ) im )dv r dv c dv s dv b (15) 

where l,m take values from 1 .. N+3 and Einstein summation convention is assumed. This represents the relative 
probability of observing any given set of these velocities in some assumed CDM universe specified by (f3, T) and the 
normalisation as. To treat any of v c , Vb, or v s as a constraint, we have to calculate the ratio of the joint probability 
to the probability of the constraint. For example, the probability distribution function of f(v r (i),i = 1..N \ v c ) is 
given by 



dv r (i) (16) 



f(v r (i),i = 1..N | v c ) = (2n)-^ 2 (^) exp -1/2 (v, ■ v m (M- 1 ) !m - 

Any other constraint is similarly easy to express. For the CMB dipole vector we use the recent determination by 
Linewaver et al. (1996) quoted in the introduction. For the bulk flow vector we use the Mark III-POTENT estimate 
at R b = 30 /i _1 Mpc(437 kms^i 40 kms" 1 , 1=306°, b=14° Dekel 1997). We have chosen not to use the Mark 
III-POTENT at large radii, which would set a more stringent constraint, because of the problems in the calibration 
of the distance indicators recently discussed by Willick and Strauss (1998) beyond 30 h^ 1 Mpc . Within a scale of 
30 h^ 1 Mpc, however, Mark III velocities have proved to be reliably calibrated and fully consistent with the 1.2 Jy 
gravity field (Willick et al. 1997 Davis, Nusser and Willick 1996). 

For the LG velocity residuals, we impose a constraint on their magnitude v s < 200 kms -1 , which means that the 
appropriate probability distribution has to be integrated over a sphere of radius v a = 200 kms -1 . The expressions for 
the various conditional distribution functions involving the first N + 2 vectors, f(v r (i) | constraints), have all been 
derived by S92, and can be easily extended to include the bulk velocity vector. 

By construction, the various probability distributions already account for errors arising from the use of finite 
window functions. However, as we have already pointed out, there are other error sources that we need to account 
for explicitly: 

• As mentioned above, we adopt the shot noise definition of S92 (their eq. (35)) which accounts for sparse sampling 
as well as the unknown mass function of the IRAS galaxies. The latter is quantified by the mass variance which for 
PSCz galaxies, in the S92 formalism, is K = 1. As displayed in figure 3, the PSCz shot noise increases with radius. 
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We calculate the shot noise contribution to the differential dipole in each shell, square it and add it to the diagonal 
elements of the covariance matrix, Ma, with i < N. 

• The PSCz dipole is measured from the real space position of galaxies which we obtain after minimising the 
redshift space distortions iteratively. Intrinsic errors in the reconstruction procedure generate an average random 
error of 15 kms _1 in the differential dipole (c/.§3.4) produced by a shell of matter 10 /i -1 Mpcwide. This error, a, is 
added in quadrature to the first i = 1...N diagonal elements of the covariance matrix just like the shot noise. 

• Nonlinear motions spoil the alignment between v r and v c . To account for these effects, we follow S92 and define 
a decoherence function in fe-space which expresses the misalignment between the Fourier components of the LG 
velocity v r ,fc and the CMB dipole v Cj fc: 

W ^ = ivn^r- < 17 > 

This additional window function multiplies the integrands that describe the covariance between measured accelerations 
and LG velocity, i.e. the terms Mij with i < N, j = N + 1 and their conjugates. Clearly, since nonlinear effects are 
more severe on small scales, Wd{k) — > at large k. To model Wd{k) we use the S92 expression 

W d (k) = 1 (18) 

(l + (fcr d ) 4 )2 

We put the decoherence length rn = 4.5 /i~ 1 Mpc which S92 templated using N-body simulation of a standard CDM 
cosmology. In adopting this expression we follow the S92 assumption that the decoherence function does not change 
appreciably with the background cosmology. 

/.From the probability density distributions introduced above and the measured PSCz dipole, it is possible to 
evaluate the relative likelihood, £, of different world models (/3, F, erg). We define the likelihood function as 

£ = -21n(/), (19) 

which is conveniently distributed like x 2 around its minimum. In the next two sections we will perform two different 
likelihood analyses with the aim of constraining the /3 and F parameters that, along with erg, specify the underlying 
cosmological model. 



4.3 Probability of Amplitude and Misalignment 

The aim of this first section is to afford us some insight into what kind of dipole we are likely to measure given a 
model of the universe and certain constraints. We do this by measuring the probability that an LG-like observer 
carrying out a PSCz-like survey measures a velocity v r in a given (/3, F, erg) universe. We characterize the LG velocity 
by its amplitude measured within a specified window, v r , and the misalignment to the CMB dipole 9. Since we want 
to plot the values of the likelihood distribution only as a function of v r and 8, we compute the constrained probability 
rather than the joint one. It is worth noting that v r scales with the assumed f3 and that the differential element in 
the probability distribution becomes dv r = ff'vl sm(8)dv r d8d(j). 

Figure 7 shows 68 %, 90 %, and 95 % likelihood contours from six probability distributions. In all cases the 
acceleration is measured within the same spherical window function of radius 150/i -1 Mpc which means that we will 
be using the window function of eqn. (^) in which R ma x = 150 h~ x Mpc . In figure 7a we plot £(v r | v c ) for a universe 
with r = 0.5, /3 = 1 and erg = 0.52, only accounting for shot noise errors (it can been shown that in the absence 
of shot noise and with an infinitely deep window function the likelihood contours reduce to a point around 8 = 
and v r = 620 kms -1 ). The variance in the predicted direction and amplitude (i.e. the width of the contour lines) 
is remarkable which shows how misleading it can be to evaluate /3 from a direct comparison between v r and v c . It 
becomes slightly larger when accounting for the errors in the reconstruction procedure, a, as shown in figure 7b. 
The inclusion of the decoherence has an even more dramatic effect on the likelihood contours. Figure 7c shows how 
seriously our ignorance in modelling nonlinear motions can affect the determination of the true dipole. The effect of 
decoherence is to eliminate the large fc-modes which greatly increases the variance in the measured dipole direction 
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(a) a B = 0.87 



(b) ct 8 = 0.52 




(c) a g = 0.87 + bulk flow 
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(d) o g = 0.52 + bulk flow 




Figure 8. Probability of (0, T). The plots 8a and 8b show the likelihood contours of £(v,-(i), i = 1, 15, v c | v s < v a ), while the 
bulk flow constraint is introduced in figures 8c and 8d. The case for a high normalisation (as = 0.87) is shown in figures a and 
c. A smaller value of as = 0.52 is used to normalise the model plotted in figures b and d. Vertical dashed lines delimitate the 
range in which 0.4 < h < 0.8 as suggested by observations. They have been drawn assuming that the normalisation is meant to 
reproduce the cluster abundance. 



and brings the most likely misalignment up to ~ 30° . The LG velocity residuals constraint imposes an upper limit on 
the amplitude of nonlinear contributions to peculiar velocities that originate on small scales and thus should reduce 
the variance introduced by the decoherence. As shown figure 7d, however, the effect on the likelihood contours is 
quite minor. Much better results are obtained when including the bulk flow constraint, as is shown in figure 7e where 
we have replaced the velocity residuals constraint with the bulk flow one. The variance in 9 is greatly reduced and the 
peak of the likelihood moves towards larger values of v r . This behaviour derives from imposing that v r and v c be well 
aligned which, together with the CMB dipole amplitude constraint, forces a significant part of the LG acceleration 
to originate within 30 /i _1 Mpc. 

Finally, in figure 7f we have changed the background cosmology (f3, T, as). The plot shows the same case as figure 
7c but refers to a as = 0.87 normalization with V = 0.195 and (3 = 0.55. This model has smaller power on small 
scales, which considerably reduces the spread of the likelihood contours. The position of the central peak, however, 
does not change appreciably. This reflects the fact that CDM models normalised to the cluster abundance produce 
large scale motions with similar properties (Jenkins et al. 1997). In particular the most likely velocity of a LG-like 
observer remains remarkably constant. 



4.4 Probability of Y) 

In the second part of this analysis, we consider a series of different (/3, V) universes and compute, for each of them, 
the likelihood of obtaining the same dipole as measured from the PSCz subsample within a specified window function 
and for a fixed as normalisation. The likelihood contours can then be displayed in a two-dimensional (/3,r) plane. 
Note that the Eke, Cole and Frenk (1996) normalisation allows one to relate as to fl m , which defines a (/?, T) — > 
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(b, h) transformation. In the following, we will, however, normally disregard this transformation and consider the 
normalisation simply as a fixed parameter. Also, allowing for a non- vanishing cosmological constant does not change 
the overall picture, as we have tested by comparing figure 8a with analogous contours produced for an Q, m = 0.3, 
Qa — 0.7 scenario. We therefore limit our discussion to the case of a 51a = universe. 

Figure 8 shows likelihood contours of /3, l eq = 1/r for two different choices of as (we have chosen to use the 
parameter l eq to allow easier comparison of our results to those of S92). Confidence levels in this case are 68, 90, 99 
%, which translates into AC of 2.3, 4.6, and 9.2 since there are 2 degrees of freedom in the distribution. The function 
plotted is £(v r (i),i = 1...JV, v c | v 3 < v a ) in which we use N = 15 non-overlapping top-hat windows of width 10 
h^ 1 Mpc . This makes use of the data we have on the growth rate of the dipole up to a distance of 150 /i -1 Mpc , as 
plotted in figure 5. 

Figure 8a,b show confidence contours of the function for both normalisations (as = 0.87 and as = 0.52 respec- 
tively). Shot noise, intrinsic reconstruction errors, a, and decoherence effects are all taken into account. Note that 
with respect to the analysis of the IRAS 1.9Jy dipole by S92 (c/.their figures 13 and 14) our constraints on F are 
somewhat weaker. This difference arises in part from the fact that, as we have pointed out, the PSCz dipole receives 
non-negligible contribution from mass inhomogeneities beyond 40 /i _1 Mpc, implying that models with large power 
(large l eq ) cannot easily be ruled out, whereas the data used by S92 does not indicate any appreciable increase beyond 
40 h^ 1 Mpc . The two different contours in figures 8a and 8b have basically the same shape, which indicates that the 
analysis is not overly sensitive to the normalisation of the power spectrum. Both exclude very small (3 (< 0.4) at 
the highest confidence level. This is not surprising since a very small /3 would not be able to produce the measured 
velocities. The range of allowed f3 is larger in the low-normalisation case and slightly shifted towards higher values 
because small fluctuations in the density contrast can only cause the observed velocities when coupled with higher 
masses. In neither case can we make any restrictive statement about Y, except that the power cannot be confined 
to the very smallest scales. The peak of the likelihood distribution is shifted towards slightly higher powers in figure 
8b, again because if the fluctuations are small and the density is generally larger, then power on large scales is more 
likely to be responsible for the observed velocities. The peak of the likelihood distribution in both cases is roughly 
consistent with the F ~ 0.2 measured by Tadros & Efstathiou (1995) from the IRAS survey. 

As mentioned above, the normalisation can be regarded as a constraint on Q. m (f2 m = 0.3 for case 8a and 
fl m = 1.0 for 8b). In that case, the /3- transforms to a b-axis and the V- to a /i-axis. In figures 8a, b the vertical lines 
indicate the range of currently accepted h ([0.4;0.8]), thereby indicating an extra constraint on F if Q. m is considered 
fixed. In both cases, this additional constraint is consistent with the likelihood contours. In case 8b, however, it only 
allows fairly high values of Y, i.e. the high fi m (coupled with a low normalisation to reproduce the observed cluster 
abundance) and a value of h in the given range force the power to originate on much smaller scales than seems most 
likely from the measured velocities (or indeed than was measured from the IRAS survey as mentioned above). 

Figures 8c and 8d show the effect of including the bulk flow constraint in figures 8a and 8b, respectively. The 
bulk velocity is mainly determined by the shape and the amplitude of the power spectrum on large scales. It is not 
surprising, therefore, that imposing the Mark III-POTENT bulk velocity constraint tends to favour cosmological 
models with excess power on large scales. Since the bulk flow is well aligned with the dipole, it also excludes the 
combination of low f3 and high F in both normalisations. In both cases, if the power is on small scales, the observed 
bulk flow can only be generated with a very high /?. However, in figures 8c,d as well as in 8a, b, the likelihood analysis 
cannot set stringent limits on the shape parameter. We only find that F < 0.4 in case of low normalisation and that 
r < 0.5 for the high-normalisation case. Both limits refer to the 68 % confidence limit level. 

The likelihood analysis allows determination of /3 with much better accuracy than F. In the high-normalisation 
case we obtain /3 = 0.71^2 with 68% confidence. Note that this is the range allowed by the likelihood contours in 
figure 8. We can tighten these constraints by marginalising over the values of F allowed by the ft-constraints. This 
yields /3 = 0.7t° 3 2 5 (see figure 9 for the marginalised probabilities corresponding to figure 8). The resulting bias 
parameter for IRAS galaxies is b — 1.43lg'J 5 , in accordance with the range of the values allowed by independent 
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Figure 9. Marginalised likelihood distribution. Marginal distribution of the likelihood of j3 obtained from likelihood distribution. 
Plots a,b,c and d refer to their analogous cases in figure 8. Continuous lines show the result of integrating over all the values of 
T. Dot dashed lines are obtained when limiting the integration range between the vertical dashed lines displayed in figure 8. 



observations and analyses. With as, = 0.52 the likelihood contours shift upwards, increasing the most likely (3 value 
to 1-01q 3 a g am at 68% confidence level. 



5 CONCLUSIONS 

We have measured the LG acceleration vector using the IRAS PSCz galaxy distribution in real space. Three in- 
dependent methods have been used to correct for redshift space distortions and to keep systematic reconstruction 
errors under control. As pointed out by S92, apart from the shot noise, the Kaiser rocket effect represents the most 
serious uncertainty in the dipole measurement. In this work we have accounted for this and all the other sources of 
systematic errors by using PSCz mock catalogues from N-body simulations. 

The resulting PSCz acceleration vector appears to receive a non-negligible contribution from scales larger that 
40 h^Mpc, unlike in the IRAS 1.9 Jy (S92) and ORS (Strauss 1996) analyses. We find that between 60 and 140 
ftT 1 Mpcthe dipole amplitude increases by ~ 35% and that only beyond this scale we find evidence for convergence. 
This result is in agreement with the recent analysis by Rowan-Robinson et al. (1998) and also confirms the reality 
of the large-scale contributions to the LG accelerations advocated by Scaramella et al. (1991), Plionis & Valdarnini 
(1991) and Branchini & Plionis (1996) on the basis of the distribution of the Abell-ACO clusters. The misalignment of 
the measured with the CMB dipole is remarkably small (~ 15°) and stays almost constant beyond 40 h^ 1 Mpc which 
corroborates the impression that most of the LG acceleration is generated from density fluctuations along the Perseus 
Pisces-Great Attractor-Shapley concentration baseline. The IRAS 1.2 Jy dipole appears to be fully consistent with 
the PSCz if the shot noise contribution is properly accounted for. The increasing importance of shot noise and rocket 
effect beyond the depth of our mock catalogues leads us to limit our PSCz dipole analysis to within 170 hT 1 Mpc. 

The measured PSCz cumulative dipole has been fed into a likelihood analysis with the aim of constraining the (3 
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parameter and the shape factor, F, assuming a CDM framework and a given normalisation for the power spectrum. 
The formalism accounts for shot noise errors, finite window size, nonlinear effects and random uncertainties introduced 
when minimising redshift space distortions. To better determine f3 and F, observational constraints on the LG velocity 
residuals, bulk flow, and LG velocity relative to the CMB frame have been introduced in the analysis. We account 
for the freedom in the spectrum normalisation by exploring two different cases: as = 0.52 and as = 0.87. In both 
scenarios we find that the PSCz dipole does not allow to set stringent limits on the shape parameter F. However, if 
the normalisation is intended to reproduce the observed cluster abundance, then a relation between as and f2 m results 
and further limits on F can be set by constraining the Hubble constant within the range allowed by observations. As 
a result we find that in a critical universe (i.e. for as = 0.52), F is constrained to be between 0.4 and 1.0 which is 
well outside the range of recent F estimates, like the one coming from power spectrum analyses of the IRAS survey 
(Tadros & Efstathiou 1995) for which T ~ 0.2. It should however be noted that the introduction of the constraints 
on h implies that T = Q m h, and will not hold for more general expressions of the shape parameter. 

The likelihood analysis produces more stringent constraints on the /3 value. For the high normalisation case 
(figure 8a) we find that f3 = 0.70^q'2 at l-a level. A result that does not change when the bulk flow constrain is 
introduced (figure 8c). These results indicate that, if cluster normalisation is assumed, then a low density universe is 
very much consistent with observations and allows determination of the f3 parameter with ~ 35% accuracy . 

When adopting a low normalisation, the constraints over f3 change appreciably with those on h and on the bulk 
flow. In particular, they change drastically with h when considering a low normalisation together with the bulk flow. 
The low normalisation implies that /? has to be very high (/3 ~ 1.5) to produce the measured bulk flow (if h is fixed 
within the given bounds). This is much higher than j3 determined from other dynamical measurements. Together 
with the fact that the range in h also forces F to be much higher than indicated by independent measurements, this 
implies that the model is probably not realistic. 

Our ability in determining (5 from the galaxy dipole is mainly hampered by our ignorance of the large scale 
contribution to the LG acceleration by the fact that the likelihood analysis is based on one single observer. Better 
determinations of /3 have been obtained by comparing Marklll POTENT measured velocities to the ones modelled 
from the 1.2 Jy catalog (Davis, Nusser and Willick, Willick et. al 1997, Willick and Strauss 1998, da Costa, Nusser 
ct al 1998). A substantial improvement, however, can only be obtained by comparing the PSCz gravity field with 
observed velocities at many independent locations. An analysis of this kind using the recent SFI data (Giovanelli et 
al. 1997a, 1997b) is currently in progress. 
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